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COMPLETE CORRECTED DIFFUSION APPROXIMATIONS FOR 
THE MAXIMUM OF A RANDOM WALK 

By Jose Blanchet and Peter Glynn 

Harvard University and Stanford University 

Consider a random walk (S n :n>0) with drift — n and So = 0. 
Assuming that the increments have exponential moments, negative 
mean, and are strongly nonlattice, we provide a complete asymptotic 
expansion (in powers of /i > 0) that corrects the diffusion approxima- 
tion of the all time maximum M = max n >o 5Vi- Our results extend 
both the first-order correction of Siegmund [Adv. in Appl. Probab. 11 
(1979) 701-719] and the full asymptotic expansion provided in the 
Gaussian case by Chang and Peres [Ann. Probab. 25 (1997) 787-802]. 
We also show that the Cramer-Lundberg constant (as a function of 
fi) admits an analytic extension throughout a neighborhood of the 
origin in the complex plane C. Finally, when the increments of the 
random walk have nonnegative mean \x, we show that the Laplace 
transform, exp(— bR(oo)), of the limiting overshoot, R(oo), can 
be analytically extended throughout a disc centered at the origin in 
C x C (jointly for both b and /i). In addition, when the distribution of 
the increments is continuous and appropriately symmetric, we show 
that E^St [where r is the first (strict) ascending ladder epoch] can 
be analytically extended to a disc centered at the origin in C, gen- 
eralizing the main result in [Ann. Probab. 25 (1997) 787-802] and 
extending a related result of Chang [Ann. Appl. Probab. 2 (1992) 
714-738]. 

1. Introduction. Let (X n : n > 1) be a sequence of independent and iden- 
tically distributed (i.i.d.) random variables (r.v.'s), and let S = (S n :n>0) 
be its associated random walk (so that So = and S n = X\ + • • • + X n for 
n > 1). This paper is a contribution to the literature on development of high 
accuracy approximations to the distribution of the maximum r.v. 

M = maxjSVj : n > 0}. 
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Clearly, —\i = EX\ must typically be negative in order that M be finite- 
valued. We shall require this condition throughout the remainder of this 
paper. The distribution of M is of importance in a number of different 
disciplines. 

For x > 0, {M > x} = {t(x) < 00} , where r(x) = inf{n > 1 : S n > x}, so 
that computing the tail of M is equivalent to computing a level crossing 
probability for the random walk S. Because of this level crossing interpre- 
tation, the tail of M is of great interest to both the sequential analysis 
and risk theory communities. In particular, in the setting of insurance risk, 
P(t(x) < 00) is the probability that an insurer will face ruin in finite time 
(when the insurer starts with initial reserve x and is subjected to i.i.d. claims 
over time); see, for example, [1]. 

The distribution of M also arises in the analysis of the single most im- 
portant model in queueing theory, namely, the single-server queue. If the 
inter-arrival and service times for successive customers are i.i.d. with a mean 
arrival rate less than the mean service rate, then W = (W n : n > 0) is a pos- 
itive recurrent Markov chain on [0,oo), where W n is the waiting time (ex- 
clusive of service) for customer n. If Woo is a random variable having the 
stationary distribution of W, then Kiefer and Wolfowitz [9] showed that 
has the distribution of M for an appropriately defined random walk. As a 
consequence, computing the distribution of M is of fundamental importance 
to queueing theorists. 

Since W is a positive recurrent Markov chain, the distribution of M can 
be computed as the solution to the equation describing the stationary dis- 
tribution of W . This linear integral equation is known as Lindley's equation 
(see [12]) and is of Wiener-Hopf type; it is challenging to solve, both analyt- 
ically and numerically. As a result, approximations are frequently employed 
instead. One important such approximation holds as fj, \ 0. This asymp- 
totic regime corresponds in risk theory to the setting in which the "safety 
loading" is small (i.e., the premium charged is close to the typical pay-out 
for claims) and in queueing theory to the "heavy traffic" setting in which 
the server is utilized close to 100% of the time. Thus, this asymptotic regime 
is of great interest from an applications standpoint. Kingman [10] showed 
that the approximation 

(1) P(M > x) w exp(-2^/cr 2 ) 

is valid as \i \ 0, where a 1 = Var(Xi). (A more precise statement of this re- 
sult will be given in Section 2.) Because the right-hand side of (1) is the exact 
value of the level crossing probability for the natural Brownian approxima- 
tion to the random walk S, (1) is often called the diffusion approximation 
to the distribution of M. 

As with any such approximation, there are applications for which (1) de- 
livers poor results. Siegmund [15] therefore proposed a so-called "corrected 
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diffusion approximation" that reffects information in the increment distribu- 
tion beyond the mean and variance. This corrected diffusion approximation 
computes the next term in the asymptotic (as fj, \ 0) beyond that given 
by the right-hand side of (1). The main result in this paper (Theorem 1) 
is a development of the full asymptotic expansion initiated by Siegmund. 
We compute all the terms in the asymptotic expansion for general random 
walks with increments having exponential moments; see Section 6 for de- 
tails on the calculation of the relevant coefficients in the expansion. Our 
theorem can be viewed as a non-Gaussian counterpart to the corresponding 
expansion provided recently by Chang and Peres [7] for Gaussian random 
walks. As perhaps expected, the mathematical approach followed here is 
quite different from that used by Chang and Peres. 

In addition to risk theory and queueing theory, the types of results and 
techniques studied here can be useful in other settings such as statistical 
sequential analysis, which is the discipline that motivated the initial devel- 
opment of Siegmund [15]; see also [16]. Other areas that have benefited from 
corrected diffusion approximations also include inventory theory; see [8] and 
computational finance; see, for example, [3]. 

As is well known in the literature, there is a close connection between such 
corrections and asymptotic expansions for the moments of the ascending lad- 
der height random variables associated with the random walk. Theorem 2 
establishes an asymptotic expansion for the mean of the first strict ascending 
ladder height for random walks with light-tailed symmetric and continuous 
increments. As indicated in Section 6, this permits one to develop asymp- 
totic expansions for all the moments of the ascending ladder heights (and 
for the limiting overshoot induced by the associated renewal process); see 
also Theorem 4. Additional expressions on boundary crossing problems and 
ladder height quantities are obtained by Lai [11] and, in a different setting 
in which fj, is held fixed, by Lotov [13]. 

This paper is organized as follows. The main results are described in 
Section 2. A key connection to asymptotic expansions for the "short-time" 
behavior of the Cauchy process is made in Section 3. Section 4 shows how 
all the integrals required for our asymptotic expansion can be reduced to 
the short-time asymptotics of Section 3. Finally, Section 5 provides rigorous 
support for the remaining details in the argument used to compute the 
coefficients in the expansion. Section 6 summarizes the computation of the 
coefficients, and discusses an expansion related to the moments of the strict 
ascending ladder height. Any proof that does not follow the statement of 
the result can be found in our final section, namely, Section 7. 

2. The main results. To state our main results, we adopt the parame- 
terization utilized by Siegmund [15]. We assume throughout this paper that 
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the Xi's have exponential moments, so that E exp(9X±) < oo for 9 in a 
neighborhood containing the origin. For such 9, define 

V>(0)=log-E(exp(0Xi)). 

Then, for each such 9, we can define the probability measure Pg having the 
property that, for n > 0, 

P e (A)=E(eM0S n -niP(9))t A ) 

for A G <r(Sj : < j < n). As is well known, S is again a random walk with 
i.i.d. increments under Pq, having common increment distribution 

P {Xi G dx) = exp(#x - ^{9))P(X 1 G dx) 

for x G R [with mean = ^'(0) and variance Vare(Xi) = With- 

out any loss of generality, assume that EX\ = and Var(Xi) = 1. Since ip(-) 
is strictly convex on its domain of finiteness, EqX\ < for < 0. Thus, Pg 
induces a random walk with negative drift when 9 < 0. We therefore focus 
on corrected approximations to Pg(M > x) as 6 /* 0. 

A key step to the analysis of Pg(M > x) is the judicious application of 
Wald's likelihood ratio identity; see, for example, [16], page 13. For 9q in 
some interval of the form (—77, 0), there exists a positive 9\ such that ip(9o) = 
ip(9i). Set A = 9\ — 9q. Note that parameterizing in terms of A is essentially 
equivalent to parameterization in terms of 9q [or parameterization in terms 
of the drift (i = —ip'(9o)]. The likelihood ratio identity then asserts that 

Pb u {t{x) < 00) = E 8l exp(-(6>i - 9 )S T ^) 

= exp(-(0i - 9 )x)E ei exp(-(0! - 9 )R(x)), 

where R(x) = S T r x j — x is the so-called "overshoot" at level x. 

Suppose now that X\ is strongly nonlattice, in the sense that, for each 
5>0, 

(3) inf |1- 9 (A)|>0, 

|A|>5 

where g(X) = E exp(i\Xi) is the characteristic function of X\ (under Pq). 
Applying renewal theory to the random walk at strictly increasing ladder 
epochs establishes then 

(4) E dl exp(-(0i - 9 )R(x)) -> E 9l exp(-(^ - 9 )R(oo)) 
as x — > 00. 

Siegmund [15] showed that the renewal theorem can be applied uniformly 
for A < 77 (see also [6]). Hence, (2) yields 



(5) Pe (M >x) = ex.p(-Ax)E 9l exp(-Aii(oo)) + o(exp(-(A + r)x)) 
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for some r > (uniformly in 8q > ~ r l/^)- I n insurance risk theory, A is called 
the "adjustment coefficient" and the quantity Eg 1 exp(— AR(oo)) is known 
as the Cramer-Lundberg constant [1]. 

Relation (5) may alternatively be written as 

(6) Pe (AM >x)=exp(-x)E dl exp(-AR(oo)) + o(exp(-rx/A)), 

where o(exp(— rx/A)) is uniform in 9q > —rj/2. Note that exp(— x) is pre- 
cisely the level crossing probability of level x/ A for a Brownian motion with 
drift —A/2 and unit variance. Since Eq x X\ ~ —A/2 as 9o /" 0, (6) provides 
rigorous support for the diffusion approximation (1). Furthermore, a cor- 
rection to the diffusion approximation described in the Introduction can be 
obtained by developing an asymptotic expansion for Eg 1 exp(— Aii(oo)). 

Siegmund [15] obtained his corrected diffusion approximation by showing 
that 



Note that, by computing the single integral (8), Siegmund's corrected diffu- 
sion approximation to the distribution of M provides a parametric approx- 
imation that is valid for all random walks having negative drift sufficiently 
close to zero. Such parametric approximations are convenient in many ap- 
plications settings (i.e., in studying the behavior of a queue when utilization 
is close to 100%). 

Our main theorem shows that there is a full asymptotic expansion for 



Theorem 1. Suppose that X\ has exponential moments and is strongly 
non-lattice. Then, r(-) (initially defined on [0,v) for some v > 0) admits an 
analytic extension on a neighborhood of the origin in the complex plane. 

Remark. An immediate consequence of Theorem 1 and the implicit 
function theorem is that the Cramer-Lundberg constant, namely, exp(r(A(#o))), 
initially defined for all 9q <0 sufficiently close to zero, admits an analytic 
extension on a disc containing the origin in the complex plane. 

According to Theorem 1, 



(7) E 01 exp(-A#(oo)) = exp(-A/3 1 ) + o(A 2 ) 



as A J. 0, where f5\ can be computed explicitly as 




r(A) = log-E^ exp(-A J R(oo)). 



(9) 
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where (5\ is given by (8) and P2 = 0. [This latter equality follows from the fact 
that the error term in (7) is o(A 2 ).] Obviously, in order for (9) to be useful 
from an applied standpoint, we need a means of numerically computing 
the Pn's. This issue is discussed in Section 6. We establish there that the 
/3 n 's can be successively computed via a finite number of one-dimensional 
integrations reminiscent of the integral appearing in (8). Thus, the /3 n 's 
can easily be computed, thereby yielding cheaply computable high-order 
parametric corrections to the diffusion approximation (1). 

The argument above also permits us to establish asymptotic expansions 
for certain ladder height quantities. As noted earlier, renewal theory applies 
to the random walk when sampled at strictly increasing ladder epochs. The 
renewal theorem invoked above actually establishes that 

, > „ / . 1 — Eo, exp(— AS T , ) 

(10) E 01 exp(-A J R(oo)) = - 



AE 9l S. 



where r + = inf {n > 1 : S n > 0} is the first (strict) increasing ladder epoch 
(see [2], page 221). In view of (2), it follows that 

(11) l-E ei exp(-AS T+ )=P eo (r + = oo). 
Random walk duality (see, e.g., page 173 of [16]) implies that 

(12) Pe (r + = 00) = 1/E 6o t-, 

where r_ = inf{n > 1 : S n < 0}. If the X^s are symmetric r.v.'s with common 
continuous distribution function, A = 28 1 and Eq t_ = Eq 1 t + . Furthermore, 
(10) to (12) then imply that 

1 

E 0l exp(-Ai?(oo)) 



28 1 (E dl S T+ )(E dl T + Y 
In view of Wald's identity, we then obtain the relation 

<//(0i) 



E 9l exp(-A J R(oo)) 



28 1 {E 01 S T+ y 



As a consequence, Theorem 1 then yields a full asymptotic expansion for 
the expected ladder height E 9l S T+ . We record this result as our Theorem 2. 

Theorem 2. Assume that X\ has exponential moments and is symmet- 
ric with a continuous distribution function. Then, 



DC 



m=0 
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Given our above argument, the only remaining issue in proving Theorem 2 
is establishing that fan — for n > 1 in the presence of symmetry. This fact 
is proven in Section 7. 

The most important device that we use to prove Theorems 1 and 2 is 
a convenient representation for r(A). This representation, is a key idea in 
our mathematical development. To introduce our representation, put (f>(6) = 
Eexp(9Xi) for 6 € K and, for zeC, set -y(z) = Eexp(zXi). Note that <j> is 
finite-valued on a neighborhood N of the origin and 7 is analytic on the 
strip {x + iy : x G Af, y 6 R}. For nonnegative 8 &Af and put 

p(6,b) = log E e exp(-bR(oc)). 

Note r(A) = p(0 x , A), where 6>i = 6>i(A) > 6> (A) = 6> is such that ^(0i(A)) = 
V>(6> (A)). Woodroofe [17] showed that 

v ' Hy ' 27ri-oo (b + i\)iX & \ -i<p'(6)\ J 

see also Corollary 8.45 and Theorem 8.51 of [16]. While (13) is convenient for 
many purposes, it presents difficulties in the current circumstances because 
of the singularity (in the logarithm) that arises when 9 \ 0. The following 
representation for p(6,b) is free of such singularities. 

Theorem 3. Suppose X\ has exponential moments and is strongly non- 
lattice. Then, for nonnegative 9 G M and b > 0, 

(14) P{ °' h) = ^ L JbTixJ-X bg { \(\-2i<W)) ) 

Siegmund's computation of f3\ takes advantage of the fact that the first- 
order behavior of r(A) should match that of 

(15) s(A) = log£ exp(-Ai?(oo)). 
Since s(A) = p(0, A), Theorem 3 implies that 

1 /-oo A 

(16) 8(A) = — / log(2(l - 5(A))A~ 2 ) dA; 

27r (A + ^A)^A 

see also page 226 of [16]. We proceed to analyze p(9,b) by writing p(6,b) = 
s(b) + 1(6, b). In view of both Theorem 3 and (16), 

(17) 7(0, 6) = - ^ log (A _ 2me)){i _ g{x)) j dA. 

In the next sections we develop asymptotics, as & \ 0, appropriate to the 
integrals arising in (16) and (17). Such asymptotics can be used to provide 
asymptotic expansions for the moments (or, equivalently, the cumulants) of 
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the limiting expected overshoot r.v. R(oo) under Pg as 9 \ 0. Specifically, 
for n > 1, let 



6=0 



Theorem 4. Assume that X\ has exponential moments and is strongly 
non-lattice. Then (for all n>\) K n {-), initially defined on [0,v) for v > 0, 
can be extended to be an analytic function throughout a disc in the complex 
plane containing the origin. 

An important implication of Theorem 4 is that it can be directly applied 
to obtain complete asymptotics for the steady-state mean of the waiting time 
sequence, namely, Eg M (= £ , g VF oo ). In particular, Siegmund [15] shows (see 
also page 270 of [2]) that 



Eg M 



(18) 



E 8q(S t+ \t + < CO) 

Po (t+ = oo) 

E ei S T+ exp(-A5 T+ ) 
l-E 9l exp(-A5 r+ ) 

E 0l (1 - fl(oo)) exp(-Ai?(oo)) 
A£ 01 exp(-Ai?(oo)) 

1 15 , . 
A + A^' A) - 



Thus, since 



-^fe)=E(-l)X + i(^, 

m=0 

it follows that Theorem 4 can be applied directly to provide the full asymp- 
totic expansion for Eq M . Indeed, our analysis in Sections 3 to 5 yield an 
asymptotic expansion for n n {-) around zero, which in turn implies the ex- 
pansion 

i n ri—m a ( AV A m 

Ee M (- 1 ) m ^ ) + i(°) M F L =T + °( An+1 ) 

Z-A n ■ n / • III. 

valid for all n > 0. The explicit computation of the derivatives 



=0 



for j, m > 0, is discussed in Section 6.2. 

Finally, the analytic extension of n n (-) and r(-) is a consequence of the 
following result. 
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Proposition 1. If X\ has exponential moments and strongly nonlattice 
distribution, then, /(•) [defined as in (17) on a domain containing [0,v) x 
[0,v) with v > 0] can be analytically extended throughout a disc containing 
the origin in C x C. 

Moreover, with the aid of Theorem 1, it follows easily [from (18) and the 
implicit function theorem] that AEg M (initially defined for 9q < 0) can be 
analytically extended [as a function of A(9q)] in a neighborhood of the origin 
in the complex plane. 

3. Short-time asymptotics for the Cauchy process. The approach de- 
scribed in Section 2 suggests computing an asymptotic expansion for r(A) 
by developing appropriate expansions for s(A) and 7(#i,A). In this section 
we will show how asymptotics for s(A) can be obtained. Section 4 shows 
how asymptotics for 1(8, A) [and, as a result, also for I(8i(A), A)] can be 
reduced to the types of integrals considered here. 

Since s(b) is real for b positive, it follows that the integral of the imaginary 
part of (16) must vanish. Hence, s(b) equals the integral of the real part of 
(16), so that 

1 f°° b 

S{P) = 2^L 7^2 Rel °g( 2 ( 1 -5(A))A- 2 )dA 

(19) 

-hIlwTWx lmlog[1 - 9iX))dX - 

Both of the above integrals take the form 
(20) 

■ 

f(Xb)d\. 



2tt J-oo 1 + A 2 

for suitably defined /. Note that if Y = (Y(t) :t>0) is a standard Cauchy 
process [so that 1^(1) is distributed as a standard Cauchy r.v.], K(t,f) can 
then be represented as 

K(t,f) = ±E(f(Y(t))\Y(0) = 0). 

Hence, representing K(t, f) as a power series in t is equivalent to the devel- 
opment of short-time asymptotics of the Cauchy process. Such asymptotics 
are also of general analytical interest, because of their relevance to Fourier 
analysis. Integrals of the type (20) are closely related to "approximate iden- 
tities of the Fejer type"; see page 31 of [5]. 

Let £ be the space of functions /:K — ► C for which £'|/(y(l))| is finite 
and for which / is infinitely differentiable at zero. For /:1R— > C, let / be 
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the symmetrization of / defined via f(x) = (f(x) + f(—x))/2. The following 
result provides our short-time asymptotic expansion for K(t,f). 

Proposition 2. Suppose f belongs to £. Then, K(-,f) is infinitely dif- 
ferentiable at the origin and 

v -l) n / 2 / (n) (0), n even, 

K {n) (0J) = { ( _ 1) (n-i)/2 n! 1 r ( T{n _ l)/2 f_)(X)dX, n odd, 



'2vr /_„ 

where, for j > 0, Tj acts on even functions in £ as 

WW ^ • 

Furthermore, the family of linear operators (T n : n > 0) is a commutative 
semigroup, so that T n+m = T n T m m,n>0. 

Remark. Note that the even derivatives of / match those of /(•)■ One 
might therefore be tempted to write the derivatives of K{-,f) in terms of 
integrals of Tjf rather than Tjf. The problem is that Tjf typically has 
a singularity at the origin, unless the odd derivatives of / at zero vanish. 
As a consequence, the integrals defining the derivative of K(-,f) may di- 
verge if they were defined directly in terms of /. To avoid this, we use the 
symmetrization /. 

Proof of Proposition 2. The fact that T n is a linear operator and 
forms a commutative semigroup is straightforward. To obtain the formula 
for the derivatives of K(-, /), at the origin, note that K(-,f) = K(-, f) where 
/ is the symmetrization of / given by /(•) = (/(•) + /(— -))/2. Furthermore, if 
/ 6 £, then / is also in £. Observe that the dominated convergence theorem 
implies that 

as t \ 0. This motivates writing 

i poo f 

K(tJ) = /(0)/2 + — J ^ (/(A) - 1(0)) dX. 

Since E\f(Y(l))\ is finite, it follows that the above integrand is uniformly 
dominated by an integrable function for |A| bounded away from zero. On the 
other hand, /(A)— /(0) = 0(A 2 ) as A — > 0, so the integrand is also uniformly 
(in t) dominated for |A| small. Hence, the dominated convergence theorem 
yields the conclusion that 

K(t, f) = /(0)/2 + — J ^ ¥ ( /(A) - 1(0)) dX + o(t) 
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as t — > 0. In fact, 

K(t, f) = /(0)/2 + - J ^ F ( /(A) - /(0)) dA 



2vr y_oo t 2 + A 2 A 2 

+ /-oo 

/(0)/2 + — / {TM^dX-eK^f). 



I : apply (21) recursively to K(-,T\f), K(-, I2/), . . . we find that K(t,f) 
satisfies 



(T J+1 /)(A)dA 



00 



+ (-i)»+i t 2 <»+ 1 )ir((,T 2( „ +1) /), 

yielding the result. □ 

With Proposition 2 in hand, our asymptotic expansion for s(A) follows 
immediately. 

4. Reducing the analysis to Cauchy process short-time asymptotics. As 

we discussed earlier in Section 2, the backbone of our asymptotic analysis for 
r(A) is given by the relation p(9, b) = s(b) + 1(9, b). In Section 3 we studied 
how to develop asymptotics for s(b). In this section we will study how to 
reduce the analysis of the remaining term 1(9, b) to that already studied in 
Section 3. Recall that 

mb) = ±- r ,~\ v log(l-v(9,X))dX, 
2-7T J-oc (b + iX)iX 



„( M)= A U 2 W) l(0)-l(0 + i\) 



where 



X-2i(p'(9)\ X 1-S(A) 

A natural strategy that one could try to pursue now is to express the loga- 
rithm as a power series in v(9,X), followed by an expansion for v as 



(22) w(M) = $>»(iA) 



n=0 



n 



! ' 



One could then apply Proposition 2 [as for (19)] to the real and imaginary 
parts in each of the resulting integrals that would appear as coefficients for 
9 n . However, the expansion (22) requires that the function v be expressible 
as a joint power series in nonnegative powers of 9 and A. Unfortunately, the 
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presence of the term (A — 2i(f>'(9)) in the denominator of v precludes the 
existence of such a joint power series. 
To avoid this difficulty, we write v as 

„(0,A)- m "- X) 



so that 



H(8,X) = l 



x-2i(j)'{ey 

2i<j/(0) _ 7 (fl) - 7 (fl + iX) 
X l-g(X) ' 



The function H(-) is well behaved because the term 2i<p'{9)/X controls the 
behavior of (7(0) — 7(6* + iA))(l — g(X))~ 1 as A \ 0. As a consequence, H (•) 
can be smoothly defined at A = via the relation H(9,0) = 1 — <p"{9). Our 
next result describes the analytic structure of H(-). 



Proposition 3. Let D v / 2 = {z E C:\z\ < 77/2} and, for (zi,z 2 ) G 
D V/2 (D V/2 UR), putH 

H{zi,z 2 ) = 1-2 -— . 

Z2 1 - 7(^2) 

Then, for every z\ E D^^, the function H(z\,-) is analytic on D v / 2 . Sim- 
ilarly, for every z 2 € -Djj/2 Ul, the function H(-,z 2 ) is analytic on D v / 2 . 
Finally, H{z\,X) can be represented as an absolutely and uniformly conver- 
gent series, for AgM and z\ € D v / 2 , namely, 

00 k 

(23) ff(zi,A) = 5> fc (iA)§ 



fc=i 



k\ 



where hk{iX) = {^ k \iX) — ^fc)/(l — 5(A)) — (2z/ifc + i/A). In particular, this 
implies that 



sup !#(«]., A) I ->0 
AeK 



as zi — > 0. 



Remark. Note that the function H{z\,z 2 ) = H(zi,z 2 ) - H(9,0) = 
H(zi,z 2 ) — 1 + 7" (21) satisfies the same properties stated for H(-) in Propo- 
sition 3 with hj,{iX) = h^{iX) + Hk+2, this follows from the analyticity of 7(-) 
and the fact that 7 "(0) = 1. Moreover, observe that completely analogous 
analytic properties apply to the function G(zi,z 2 ) = {z\))~ l H{z\, Z2) de- 
fined on D v i 2 x (-D^/2 U R) . 
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Note that |A/(A - 2icf>'(9))\ = |A|(A 2 + (2^(^)) 2 )~ 1 / 2 < 1. It follows from 
Proposition 3 that, for r > small enough, 

sup sup \v(9, A) | < 1. 

6»G(0,r) AeK 

Therefore, for all < 6 < r, we can proceed to expand log(l — v) in powers 
of v and formally integrate each term in the obtained expansion to express 
I(9,b) in terms of integrals of the form 

(24) A(a , 6 , /) = _/_^_^_) m)dK 



(b + \a + iX, 

where a, b > 0, f(i-) G C and k > 0. Because Jo(a,b,f) = Jo{b, /) can be 
written as 

/oo z. 



Mb,f)= 1 

(25) 



it follows that asymptotics for Jo can be computed in terms of asymptotics 
for the if- type integrals that are subject of Proposition 2. In view of the 
development leading to (24), a key to our asymptotic expansion for I(6,b) 
is therefore the reduction of integrals Jfc(a, b, f) for k > 1 to integrals of 
the form Jo (6,/). A key identity in establishing this reduction step is the 
following. 

Lemma 1. Suppose that a, b > 0. Then, for m, n > 0, 

j_ r ~ b (iA)m+1 dx=o 

2vr 7_oo (b + iX)iX (a + iX) m + n + l 

Furthermore, 

1 Z" 00 -1 

T- / 7T— avY lo S 1 + oiA ) dX = °- 

27T J_oo (1 + «A)zA 

Proof. For a,b > 0, let the function of a complex variable /(•) be de- 
fined as 

f(z) _ -b (^r +1 

1 ; (fr + iz)iz (a + iz) m+n+1 ' 

Consider the contour (in the clockwise direction) C(r) = C\{r) + C2(r), 
where Ci(r) = {re"~ : — tt < t < 0} and C*2(r) = {A: A 6 [— r, r]}. Since / is 
(complex) analytic on Im(z) < 0, Cauchy's theorem yields 

1 r -" < i2 > m+1 *-<>. 



27T J C (r) (b + iz)iz (a + iz) m+n+1 
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This, in turn, implies that 

1 f r -b (iX) m+1 



2vr J- r (b + [a + iA) m+n+1 

-1 r -b (iz) m+1 



dX 

n 

dz 



27r 7ci(r) (b + iz)iz (a + iz) m+n+1 
2vr J — 7r (b + ire Ti )(a + i r e Ti ) m+n+l 



dr. 



Letting r — ► oo, we obtain (by virtue of dominated convergence) the first 
part of the lemma. For the second part, let us define 



l r°° 

/i(o) = 5-/ -((l+tA^-MogCl + aiA) 



dA. 



A routine dominated convergence argument, combined with our previous 
analysis, shows that 

i poo _1 

/'(a) = — / dX = 0. 

Jn 1 27ri_oo (l + a)(l + azA) 

The proof of the lemma is complete by observing that f\{a) — > as a \ 0. 
□ 

Let £o be the subspace of £ (recall the definition of £ preceding Propo- 
sition 2) for which /(0) = 0. Also, for / 6 £, let /(■) = /(•) - /(0) (e -Co). 
We are now ready to offer a proposition that reduces the evaluation of the 
integrals Jk(a,b,f) for k > 1 to that of integrals such as Jo(b,f), thereby 
permitting the application of the short-time asymptotics of Section 3. 

PROPOSITION 4. Suppose that f € Co- Then, for k>l and n>0, 
(26) J k (a, bJ) = J (b,J2( k + J .- 1 ^ (-ayfj^j + bo(a n ), 

where the linear operator Tj (j > 0) acts on functions f(i-) € Co as 

m)-J2L=J {m) (o)(iX) m /m\ 



(T 3 f)(iX) 



(ixy 



Moreover, the family of operators (Tj : j > 0) constitutes a commutative semi- 
group, so that f m f n = f m+n . 

Remark. As for Proposition 2, one might be tempted to express the 
right-hand side of (26) in terms of / rather that /. However, Tjf is generally 
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nonintegrable with respect to the kernel that defines Jo- Finally, note that, 
if all integrals are interpreted in terms of Cauchy principal value, one can 
apply Proposition 4 directly to functions that do not vanish at the origin by 
defining J (6, /) = J (b, /(•) - /(0)) + /(0)/2. 

Proof of Proposition 4. That (fj : j > 0) is a family of linear oper- 
ators forming a commutative semigroup is immediate. By virtue of Lemma 
1, it follows that 

J m (a,b,f) = ±- [°° ~\ (J±-Yf(i\)d\. 
27r J_oo (b + iX)iX \a + iXJ 

Observe that /(«■) is now in the domain of the operators T n , n > 1. On the 
other hand, we can write 

(27) + - r k "» jw ( (-^t) ra - i v 

lit Loo (b + iX)iX J v \\a + iXJ J 



Note that 



iX \ m 1 _ g a fc (U) m " 



a + iXJ fr[\ k J (a + «A) m ' 

Once again, by appealing to Lemma 1 and to the definition of T^f, it follows 
that, for m > k > 1, 

fe ~ j\ a 1 f°° ~b a k {iX) m " k 



2ir J-oo (b + iX)iX [a + iX) m 
Combining this observation with (27), we obtain 

J m (a,b,f) = J m (a,b,f) 

(28) 

= J (bJ)-J2 [™) a k J m (a,b,T k f). 

The recursive relation (28) can now be expressed in operator form as 

J m (a, b, J) = J (b, J) + J m (a, b, (1 - (1 + af) m )J). 

(Here, we have used the semigroup property of the family of operators T m .) 
Iterating the previous expression, we arrive at 

J m (a,bJ) = J m (a,b,f) 
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(29) 



£jo(&,(l-(l + aTr/7) 
k=0 



+ J m (a,b,(l-(l + aT) m ) n+1 f) 

-ayfjf] + bo(a n ), 



\ j=0 



m + j — 1 
j 



where the last equality in (29) has been obtained by using the semigroup 
property of the operators T m and by noting that the coefficient of a^Tj in 
(29) (for j < n) must match that of x J in the formal expansion of 



p(x) 



m\n+l 



l-(l-(l + x) m ) 
1 - (1 - (1 + x) m ) 
1 



(l + x) r 



+ 0(x 



n+l^ 



That the error term in (29) is bo(a n ) comes from the fact that aJ m (a, b, f) 
bo(l), as a \ 0, as it can be seen as follows: 



\aJ m (a,b, f)\ 





a 






2^ 


J — oc 




b 


roo 


< 








2tt J 


— oo 



-b(i\) m ~ l f(i\a) 
{b + i\a)(l+i\) m 

f(i\a) 



dX 



X(l + iX) 



d\ = bo(l), 



where the last step follows by a dominated convergence argument. This 
concludes the proof of Proposition 2. □ 



Proposition 4, combined with our development for K(t, ■) in Section 3, 
provides all the elements required to develop asymptotic expansions for in- 
tegrals of the form J m (a,b,f). Since, as discussed earlier at a formal level, 
1(9, b) can be expressed as a sum of terms such as J m (o 5 b, /), it follows that 
the whole asymptotic analysis of r(A) and p(0,b) can be reduced to that of 
Section 3. A complete rigorous justification for this representation for 1(9, b) 
is one of the main issues discussed in Section 5. 



5. An asymptotic expansion for l(0,b). In Sections 3 and 4 we have 
developed the tools required to obtain asymptotic expansions, in powers of 
b, for s(b) and 1(9, b). We have done this by showing that the problem can 
be reduced to short-time asymptotics for the Cauchy process. The purpose 
of this section is to make rigorous the expansion for 1(9, b), in powers of 9, 
that was outlined in Section 4. 
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Noting the important role that functions vanishing at the origin plays in 
Proposition 4, it seems appropriate to define 

H(9, A) = H(9, A) - H(9, 0) = H(9, A) - 1 + <t>"{6) 

= ^/i fc (iA) — , 
fc=i K - 

where h k (i\) 4 ( 7 ( fc )(a) - - 5 (A)) - (2^ fc+1 /A) + /i fc+2 is such that 

/ifc(0) = 0. The next proposition shows how a simplified expression for 1(9, b) 
in terms of H can be obtained. 

Propositions. Define -$(9) = 20 (6) /ft' (9). Then, 

(31) ^' 6) = ^ loo (6TO)a l0g I 1 (A-itt(6>)) ) dX - 

Proof. Just note that 

\H(9,X) \(l-<p"(9)) 



log(l-«(e,A)) = log(l 
= log 



(A — i^(9)4>"(9)) (\-i^(9)ft'(9)) 
»A/¥(0) + l 



i\(j)"(9)/^(9) + 1 

( f)"(9)- 1 XH(9,X) 



log 1 



(A -**(*)) 

Thus, (31) follows from Lemma 1 by noting that 



2vr y_oo (6 + t A)iA \i\<//'(9)/V (0) + 1 



2W_oc (l + «A)iA & V*A&<£"(0)/tf(0) + l 
□ 

Additional simplifications reduce the complexity of the expansion for 
1(9, b). In particular, the expression for the integral Jo(b,f) simplifies when 
it is known that Jo(b,f) is real; see (25). Fortunately, our analysis of 1(9, b) 
gives rise to such real- valued Jo(b, /)'s. To establish this result, we introduce 
the following family of functions. 

Definition. A function / : K — ► C is said to have the "parity property" 
if Re/(i-) and Im/(i-) are even and odd functions, respectively. The class 
of functions possessing the parity property will be denoted by V ■ 
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Note that if f(i-) is in the domain of Jo(b, •) and / possesses the parity 
property, then we must have that Im Jo(b, f) = (since it corresponds to 
an integral on the real line of an odd integrable function). The family of 
functions enjoying the parity property has certain closure characteristics 
that will be useful for the rest of our development. These closure properties 
are discussed in the next proposition. 

Proposition 6. The class V of functions forms an algebra on R (i.e., 
a vector space on R that is closed under product of functions). In addition, 
if f & V, then l//(-) (defined on its domain of finiteness) also possesses 
the parity property. Finally, if f is in the domain of T and has the parity 
property, then Tf 6 V . 

Proof. Certainly V constitutes a vector space on R and it is almost 
immediate that T preserves the parity property. Now, if /i,/2 € V, then 
Re(/i/2) = Re(/i) Re(/2) — Im(/i) Im(/2) must clearly be even. Similarly, 
Im(/i/2) must be odd, which implies that /1/2 € V. Finally, note that 

I = Re(/) _ Im(/) 

/ Re(/) 2 + Im(/)2 *Re(/) 2 + Im(/)2' 

which immediately implies that Rel// and Iml// are even and odd func- 
tions respectively and, thus, 1// £ V- □ 

We now present the main result of this section, which yields an expansion 
for 1(9, b) in powers of 9 and coefficients involving only integrals of the form 
Jo(b,f) with / satisfying the parity property. 

Proposition 7. For k,m> 1, let the coefficient multiplying 6 k in the 
power series representation of G(9,X) m = ((p"(9)~ 1 H(9,X)) m be defined as 
9k,m(i^)- Then, gk,m(-) £ V can be recursively computed via 

k 

<?ifc,m+i(«A) = ^2g n+ i )m (i\)g k -n,i(i>»)- 

n=0 

Consider 6 > and let x(9) = -^{9)/9. Then, 

n m—1 

(32) 1(0, b)=Y,G m T. x(9) J Mb, E hm ) + bo{9 n ), 

m=l j=0 

where Ej jTn = Ej tm (iX), defined for < j < m — 1 and m > 1 as 

TP - _ m V s 1 1 ( m - k - l\ ~ ^ 

-~3,m — / , ■ , - I 1 j9k,m— j — ki 

^ m-j-k\ J J 
satisfies the parity property. 
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Proof. Since 7(iA) = E Q cos(\Xi) + iE Q sm(\Xi), it follows that 
fj,k € V, as does the function 1 — 7(iA). By the closure properties described 
in Proposition 6, we may easily conclude that g^ i £ V . A second application 
of Proposition 6 shows that gk, m £ and Ej m € V . The recursive expres- 
sion provided for gk,m follows from standard convolution operations of power 
series. For n > 1, define 

n ak 



and 



Note that 



a 1 f°° ~b . A G n (0,A)A 



W 6) = 2^ (6TO)zA l0g I 1 " A-aV(g) J 



X-i2<t>'(9)J b V X-i2<j)'{9), 

A (G n (0,\)-G(9,\)) 



log 1 



A - »2^(fl) (1 - G n (0, A)A(A - i2^(fl))- 1 ; 

On the other hand, from the remark following Proposition 3 and because 
log(l + z) = z(l + e(z)) for zGC, where < \z\ for |zj < 1/2 (see Propo- 
sition 8.46, [4]), we can see that there exists a constant B > such that 

£. /■«> b\G n (9,X)-G(9,X)\ 

~ 2TT Loo (&2 + A 2 )V2(A2 + (2^(^)2)1/2 • 

Essentially, by making the change of variables u = X9, we then see that, for 
all 9 € (0, 5) for some 5 > 0, we have 



2vr 7_oo ^+ 1 |A|(A 2 + 1)V2 

It follows easily from the previous inequality and the dominated convergence 
theorem that 

I(9,b)-I n (9,b) = bo(9 n ). 

Using the expansion of log(l + z) at z = and a similar dominated conver- 
gence argument, we can write 

_1 roc _h n 1 / ,-\ \ m _ 

In(9,b) = — ° X G n (O,\rd0 

2tt J_oo (o + zAjzA ^ 1 m \«A + 

+ 6o(r ) 
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-1 
2^ 



E 



i 



iX 



3 (6 + iA)iA \iX + ^(0) 

+ bo(8 n ). 

Using Proposition 4 and (33), we obtain that 
1(9, b) 



m n—m 



Y.° k+m 9kA^)dX 

k=0 



n m—1 , 

-Y, em J2 J ^-kh(o),b, 

m=l k=0 ^ 

+ bo(0 n ) 

n m—1 / n 



9k,: 



m — k 



(34) 



m=l k=0 

+ bo(9 n ) 

n m—1 



j=0 



m — k + j — 1 
3 



9k, 



m—1 



m=l i=0 

+ 6o(0 n ) 

ra m—1 



fc=0 



m-i-1 



m — k + j — 1 
J 



J m — k 



9k,m 

m — k 



k=0 



m — k — 1 



9k,m(iX) 
m — j — k 



m=l j=0 \ 

+ bo(6 n ), 

which yields the desired conclusion. □ 

In view of the previous result, an explicit expression for the coefficients 
in the expansion for Jo(-,f), when / satisfies the parity property, deserves 
special attention. Providing such explicit expressions is the aim of the next 
proposition. 

Proposition 8. Suppose that /(£•) £ Co has the parity property. Then, 
Jo('if) is infinitely differentiate at zero and 



f(-ir /2 te(o) 



2^ 



I poo 



(T n /2+ifm){X) dX 



n even, 



(_i\(n+l)/2 ffiu J (0) 
1 ' I (n + 1) 

~2^J ( T (n+i)/2fRE)(X)dX), n odd, 



(35) 
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where /im(*A) = Im/(iA)A 1 and /re(«A) = Re/(iA). 

Proof. The proof follows by a direct application of Proposition 2 com- 
bined with the fact that Re Jo(b, f) = 0. □ 

We close this section with some remarks that clarify how the expansion 
just derived for 1(6, b) can alternatively be viewed through the prism of a 
formal operator expansion. The analytic properties stated in Proposition 1 
provide rigorous justification for the expansions outlined next. First, we note 
that if 9 > is small enough and b > 0, we can formally write 

(36) 1(9, b) = -J2 j;Jk(*(9),b, <f>" \e)~ k H k (6 , ■)). 

k=i K 

Formally interpreting (1 + aT)~ m as 

(l + af)- m = Y j ( m+ k~ 1 ) (~a) k f k , 
k=o ^ ' 

in combination with the expansion (26) developed for Jf : (a,b,f) and equal- 
ity (36), allows us to write 

1(6, b) = -J2 rMb, 4>"(9Y k (l + *{0)f )- k H k (9, •)). 

If we introduce the convention that, for commutative operators B\(6), 
B2(6) and functions F\(9, •), ^(0, •), expressions of the form B\(9)F\(9, •) x 
B2(6)F2(6 , ■) (or any permutation of this form) are always interpreted as 

(B 1 (9)B 2 (9))(F 1 (6,-)F 2 (6,-)), 

then we can write 

(37) 1(6, b) = J (b, log(l - <tJ\0)-\l + ^(6)fy 1 H(6, •))). 

Expression (37) provides a convenient shorthand notation for the expansion 
of 1(6, b), in powers of 9 and with coefficients in terms of integrals of the 
form Jo(b, •). In addition, note that, in order to recover the coefficients in 
the expansion for I(9i(-),-), one can apply formal differentiation to (37) in 
both arguments 9 and b [always having in mind that (37) is just a formalism 
representing a certain asymptotic expansion]. Hence, for example, one can 
obtain the first term in the expansion for I(6i(-), •) as 



d A /(0i(A), A) |a=o = d e I(0, 0)<9 A 0i(O) + d b I(0, 0), 
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where the formal derivatives applied to (37) must be interpreted using the 
formal operator convention introduced earlier. Thus, for example, if B{9) is 
an operator of the form 

oo ry<k 

k=o K 

applied to a function F(9, X) = J2 fk(iX)9 k /kl, we interpret the formal deriva- 
tive dfllog(l -B(9)F(9,-)) as 

do log(l - B(9)F(9, •)) = -d e B(9)(l - B(9)F(9, ■)y 1 F(9, •) 

-B{9){1- B(9)F(9, ■)y 1 d e F(9,-), 

where 

d e B(9)(l-B(9)F(9,-)r 1 F(9,-) 

oo 

= Y J {d e B{9)B{9) k )F{9,-) k+1 

k=0 

and, similarly, 

B(9)(l-B(9)F(9,-)r 1 d e F(9,-) 

oo 

= Y J B{9) k+ \F{9 ) -) k d e F{9,-)). 

k=0 

Thus, it is possible to combine this formalism with the expansion 

m 

J (b, f) = J2 J {n) (0, f)b n /n\ + 0(b m+1 ) 

n=l 

to recover the coefficients in the expansion for I(9i(-), •) in powers of A. 

6. Expansions for r(A) and EoR k (oo). In previous sections we de- 
veloped all the elements required to rigorously compute a full asymptotic 
expansion for r(-) in powers of A. In the first part of this section, as a 
summary, we indicate how the developments obtained in the previous three 
sections can be applied to provide an asymptotic expansion for r(-) in powers 
of A. In view of the level of complexity in the computation of the constants 
(3 n , the description in this section is intended to provide guidance for an 
easy-to-design practical implementation in a computational package such as 
Mathematica or Matlab. An efficient implementation of the procedure will 
appear elsewhere. In the second part of this section, also as a direct conse- 
quence of the analysis in the previous sections, we will develop a rigorous 
asymptotic expansion for the cumulants of R(oo) under Pq in powers of 9. 
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6.1. The expansion for r(A). An algorithm for computing (3t for k <n 
proceeds as follows: 

1. Expand s(A) up to terms of order 0(A n+1 ) using Proposition 8. 

2. Similarly, expand the functions Jo(-, Ej,m) up to terms 0(A n ~ m ) with < 
j < m — 1 and 1 < m < n. This also can be done by applying Proposition 
8, since Ej >m has the parity property. 

3. Finally, the terms obtained can be combined with an expansion for #i(A) 
up to terms of order 0(A n+1 ). Such an expansion can be easily obtained 
using the implicit function theorem and therefore is omitted. 

Observe that the previous algorithm provides an asymptotic expansion 
for r(-) in powers of A. However, because of Theorem 1, we actually have 
that this asymptotic expansion converges absolutely in a neighborhood of 
the origin. 

As a simple application of the previous expansion, we show that @2 = 0. 

Proposition 9. Suppose that X\ has exponential moments and is strongly 
nonlattice. Then 

r(A) = -A/?i + 0(A 3 ). 

Proof. We only need to show that /?2 = 0. Note that, by virtue of 
Proposition 8, the coefficient multiplying A 2 in the expansion of s(A) equals 

1 f°° 1 /Imlogfl -5(A)) 
S2 = ^jLaH A ^ 

t± _ M 

12 18 

In order to show that = 0, it suffices to show that 6\ J(A,Eq^) ~ — A 2 S2 
or [since J(A,E jtm ) = O(A), 6 1 /2 ~ A and <£"(0i) ~ 1] that AJ(A,%) ~ 
— 2A 2 S2, where 

= L wrm (ra " a + dX 
a* 1 r A2 p ( 2 w \^ 

(38) = nj WT^ Re {T^giX)-J + ^) dX 

A f°° A 2 / 2i \ 

< 39 » -7/o WTwMTrjT)-J + ^) dx - 

Note that g'(X) = i^'(iX) and that Imlog(2A -2 ) = 0; hence, we can write 



Re( i _ g (A) " 7 + ^ 3 J = " Im ^( 1 °g( 2 ( 1 - 5(A))A~ 2 ) - fi 3 iX), 
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which implies, using integration by parts, that the integral in (38) equals 

— 1 f°° 2AA 2 

— J (A 2 + A 2 )2 Im(log(2(l - g(A))A~ 2 ) - » 3 i\) dX 

(40) 



-— f 2/Imlog(l- g( A))_ x 
7r Jo A^ V A / 



where (40) has been obtained using dominated convergence and simple ma- 
nipulations. It follows from Proposition 2 and a first-order asymptotic ex- 
pansion of Eo^iX) that (39) equals — A 2 (/i 2 /9 — ^4/6). Combining this last 
estimate together with (40) into (38) and (39) yields AJ(A,%) ~ -2A 2 s 2 , 
which is exactly what we wanted to show to conclude that fa = 0. □ 

6.2. The expansion for EQR(oc) k as 9 \ 0. We shall provide asymp- 
totics for E e R(oo) k = E e (S* + ) / '(k\E e (S T+ )) via the cumulants («j(6>) :j > k) 
of R(oo) under Pg. In particular, these estimates yield the proof of Theorem 
4 stated in Section 2. The idea is to develop an asymptotic expansion, in 
powers of b, for s(b) and 1(9, b) respectively and to match coefficients in the 
expression 

p(9, b) = - Kl (9)b + K 2 (9)b 2 /2 - K 3 (9)b 3 /3\ + ■■■ 

(41) 

= s(b) + I(9,b). 

In order to perform this task, we will take advantage of Proposition 7 as 
follows; first let us define, for k> 1, ctk,j, m = Jq (0, Ej^ m )/k\ (which can be 
explicitly computed via Proposition 8). With this notation, we can write, 
for l,n> 1, 

n m—l / I \ 

1(9, b)=J2 0™ E X(0) j E a k, hm b k + 0(b l+1 ) + bo(9 n ) 

m=l j=0 \k=l I 

I n m—l 

= E bk E E m x(eya kJ , m + 90(b l+1 ) + bo(9 n ). 

k=l m=l j=0 

Therefore, we obtain that, for all s,n > 1, k s (9) satisfies 

(n m—l \ 
Ks(0)+s\ E E m x(eya sJ)m )+O(9 n+1 ). 
m=l j=l ) 

Consequently, n n (-) is an infinitely differentiable function at 9 = and for 
m > and re > 1, we have 

4" } (o) , n n.^(o) T; lm ^ s 

s=0 j=0 
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where, for n,j>l, Xn,j is the coefficient multiplying 8 n in the expansion for 
x(d) 3 . In particular, the Xn,j can be computed recursively as 

Xn, ]+ i = j2xn, 3 X (k - n) (0)/(k-n)\, 

n=0 

with Xn,i = X (n) (0)/n!. 
7. Technical proofs. 

PROOF of Theorem 3. Using Lemma 1, we can add 
1 r°° b 

to expression (13) for p(6,b) to obtain 



2vr 7_oo (6 + »A)iA V -#'(0)A(1 + i\/2(j)'(e)) 

i r -6 log p( 7 m- 7 (»+a))N^ 



2vr 7_oo (6 + V A(A-2^(0)) 

yielding the conclusion of the theorem. □ 

Proof of Proposition 3. It follows immediately, by a Taylor series 
expansion of 7(-), that a series representation for H can be written (for fixed 
A and 9 such that < |A| + \9\ < r/) as 

2icj)'(0) 7 (0)- 7 (0 + iA) 



H(6,\) = l 



A 1- 5 (A) 



fc=l 

In fact, the functions h k {i-) can be analytically extended throughout the 
disc D v j2 = {z £ C: \z\ < r]/2}. This is easily seen as follows, recall that 
7(-) [and therefore 7^(-)] are analytic on AT (defined in Section 2). Also, 
observe that 1 — j(iz) ~ z 2 /2 and j( k \iz) — fi k ~ izp k+ i as z — > 0. Thus, 
(-y(*0 (^) — — 7(iz)) possesses a simple pole at with residue equal to 

2ijLtfc_|_i , which implies that the natural extension of % defined as 



h k {iz) 



-i (k \iz)-p k 2ifi 



k+l 



1 — "f(iz) z 

(<y( k \iz) - [i k )z - 2i/i fc+ i(l - j(iz)) 
(l-j(iz))z 
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is analytic on D v /2 - Now, by virtue of the maximum principle (see, e.g., [14], 
page 253), we have that if 5 > is suitably small, 

sup \h k (iz) \ < sup \h k (iz)\. 
\z\<8 \z\=S 

Since j(z) is a nonconstant analytic function defined on £^/ 2 (which is an 
open set and thus has an accumulation point), then 1 — 7(2:) has an isolated 
zero at z = 0. Thus, it is possible to choose 5 > in such a way that 

inf |1 - y(iz)\ > e > 0, 
|*|=5 

for some e > 0. Consequently, 

sup \h k (iz) \ < sup \h k (iz)\ 

\z\<8 \z\=S 



< ~7 sup\(j^(iz)-fi k )z + 2 f i k+1 {l-- f (iz))\. 



e $ \z\=5 

Observe that, for \z\ <r//2, j( k \z) = E (X k exp(zX)). Therefore, if z = 
x + iy, with \z\ = 5, 

\^ k \iz)\<E (\X\ k \eMizX)\) 

= E (\X\ k \e W (yX)\) 
<E (\X\ k exp(S\X\)). 

A similar bound can be obtained for 7(2) and we can conclude that 3B > 
such that 

sup \h k (iz)\ < B(E (\X\ k (eMS\X\) + 1)) + E (\X\ k+1 )(l + E exp(5\X\))). 

\z\<8> 

Now, suppose that 5 < r//2. Then, if z\ G D^n, we can define 

BE (\X\ k+1 )(l + E eM$\X\))ji = Ni(zi) 

in such a way that the previous series converges absolutely and uniformly 
on -Djj/2- Similarly, we can define 

00 k / 00 k \ 

B J2 E (\X\ k (exp(5\X\) + 1))| = BE, ^ |X| fc | (exp(5 | X |) + i) 

k=l \fe=l / 

= BE ({exp( Zl \X\) - l)(exp(6\X\) + 1)) 
= N 2 ( Zl ). 
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Note that, for j = 1,2, Nj(zi) — > as z\ — * 0. On the other hand, since g(X) 
is strongly nonlattice, we have that 



sup \hk(iX) \ = sup 
|A|>5 |A|><5 



7 (fc) (i A) - nu 2i/x fc+ i 



1-5(A) A 
<B(E (\X\ k ) + E (\X\ k+1 )), 

if i? < oo is big enough. The previous estimates imply that there exist con- 
stants < M k < B{E (\X\ k {exp{5\X\) + l)) + E {\X\ k+1 )(l + E exp(5\X\))) 
such that 

\h k {iz2) \ < M k 

for z 2 € RUZya and | Y%=1 M k ^\ < EfeLi |M fc $| < oc for zi G A, /2 . Thus, 
using the Weierstrass M test, we obtain the validity of (23). Finally, the 
invoked Weierstrass M test, combined with the analytic functions conver- 
gence theorem (see Theorem 10.28, page 214, of [14]), yields the analyticity 
of H(zi, •) on (f° r z i £ ^77/2) an d similarly for H(-,Z2) on (for 

z 2 elufl, /2 ). □ 

Proof of Proposition 1. We start by writing 

The strategy will be to study this integral on {|A| < 5} and {|A| > 5} sep- 
arately (where 5 > is some convenient small number to be characterized 
later) : 

(43) - T, J lM > s (MAja log I 1 - x-2<p W ) dX - 

Let us define Ia(0,&) and Ib(0,6) as (42) and (43) respectively. Suppose 
that < 6 < 5 < 77/2. By making u = 6A, we can write 

6) = - - y_ a - log ^ - Y ^ ww ) dX. 

Let C = {w £C:\w\ < 5} D {Im(u>) < 0}, and observe that, by virtue of 
Proposition 3, we can pick S\ > in such a way that, for all < 9 < 5±, the 
function 

[b + iw)iw \ w — 2(p'(6)i/ 
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is analytic on C. Thus, applying Cauchy's theorem to the contour enclosing 
C, we obtain 



2vr J- n (b + iSe iX ) iSe iX 

H(9,8e iX )8e iX 



(44) 

1 f iW _1 e -<A 



2vr J- n (l-ib5- l e- iX ) 

Equality (44) has been obtained by simple algebraic manipulations. Observe 
that the previous expression in combination with Proposition 3 and the 
analyticity of the functions <f>'(9) (~ 0) at zero immediately gives that Ia(0, b) 
can be represented as an absolutely convergent double power series in 9 and 
b on the set < \9\ + |6| < 82 for some 82 > 0. Indeed, if we pick 82 small 
enough, it is possible to provide an explicit power series representation for 
Ia{9-, b) by using the expansion of log(l — w) at w = in combination with 
the series representation (23) for the function H (9, A) derived in Proposition 
3 and a Taylor expansion of (1 — w) -1 around w = 0. 
The analysis of lB(9,b) is easier, 



1 , ft H(9,X) 



l-i2ft(6)\- 



dX. 



Hence, in order to show that Ib{~) can be written as an absolutely conver- 
gence double power series in a neighborhood of the origin, it suffices to show 
(by Fubini's theorem) that 

~ f m + k\ V +1 (2(E(exp(9\X\) - 1 - \X\))) m 
M*kh*\ k > (fc + l)|A^ 

is finite for all nonnegative 9 and b such that 9 + b < 83 for some 83 > 0. 
But this fact follows easily from Proposition 3, first note, by the change of 
variables X = u8, that the previous expression equals 

+ k\ V +1 {2{E(exp(9\X\) - 1 - \X\))) m 



T m 
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/ oo n 

x £ \h s (i\6)\- 

\s=l S - 

now pick 63 small enough so that < max(6, 2(E(exp(9\X\) — 1 — |X|))) < 
63 < 5 (if 9 + 6 < £3), and use Proposition 3 to conclude that one £3 can be 
chosen so that \h s (i^8)\jr < c<l — 63/6. Therefore, we can bound the 

previous sum by 



2 1 

< 



3 1 - 83/5 



log 1 



1 - S3/S 



< 00. 



The conclusions obtained for both Ia(0 and -/#(•) indicate that, for all < 
0,b<v (for some v > 0), 7(0, 6) can be written as 

7(0,6) = ^ 

j,fc>i 

where the previous series converges absolutely on the specified region on 9 
and b. The previous expression provides the natural analytic extension of 
!(•) onDl = {(z 1 ,z 2 )eCxC:\z 1 \ + \z 2 \<v}. □ 

Proof of Theorem 1. Since 

exp(s(A)) = AE (S ) = £; o(exp(-A J R(oo))), 

the analytic extension of the term s(A) follows from that of the right-hand 
side, which comes from the fact that S T+ has exponential moments (see [2]). 
Thus, since r(A) = s(A) + /(0i(A), A), we just have to analyze I(0i(A), A). 
However, from the implicit function theorem, we know that 0i(O is analytic 
in the neighborhood of the origin, thus, the analytic functions convergence 
theorem (see Theorem 10.28, page 214, of [14]) combined with Proposition 
1 yields the desired conclusion. □ 

Proof of Theorem 4. From Theorem 1, we know that, for < 9, b < v 
(for some v > 0), 

00 

I{6,b) = Y,VL j {0), 

where each function I.j(9) can be expanded in absolutely convergent power 
series for < 9 < v , and thus can be analytically extended throughout a 
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neighborhood of the origin in the complex plane. But, 

p(0,b) = -ki(0)6 + K 2 {6)b 2 /2 - K 3 (6)b 3 /3\ + • • • 
= s(b) + I(6,b), 

where s(-) is (real) analytic at zero. Hence, the conclusion of the theorem 
follows immediately by matching coefficients. □ 

Next, we show that if the distribution of X± is symmetric, then, for n > 1, 

&2n = 0. 

PROOF of Theorem 2. As we discussed before, all that we need to 
show is that fan = 0. We have shown that an absolutely convergent power 
series representation is possible for r(A) when A is small, thus, it suffices to 
show that if < A < 5 (where 5 > is suitably small) , then an asymptotic 
expansion for r(A) is given in odd powers of A only. Using the integral 
expression (14), integrating on |A| < 5 and |A| > S, we can write 

(45) r{A) = 2^J lxl< s^TiX^ l0g { A(A-2 Wl )) J"" 
(An 1 f -A j 2 ( 1 (0 1 )- 1 (0 1 + l \)) \ 

(46) + J W > S JK+iXjiX bg { A(A-2 Wl) ) J ^ 

Define by ^4(A) and -B(A) the integrals appearing in expressions (45) and 

(46) , respectively. We first analyze ^4(A). Using a similar argument as in the 
proof of Theorem 1 , we see that 

A(A) — — ( — \oJ 2{l{6l) - l(6l + lZ)) \dz 
A[/X) -27rJ Cl (A + iz)iz [ ° g { z{z-2i^{e l )) ) dZ > 

where the trajectory C\ is defined as C\ = {5e lX : A G [0, — 7r)}. Also, define 
the trajectory C2 = {5e lX : A G [— n, 0)}. The proof of the theorem will be 
complete if we show that ^4(A) is an odd function. That is, we must show 
that A(A) = -A(-A). Note that 

a(-a^- zl f ~ A 1 / 2( 7 K)- 7 H 1 + K )) \ 

1 ' ~ 2tt J Cl (- A + iz)iz g \ z(z - 2i0(-e 1 )) J 

(47) 1 f -A . log p(7^) + 7(^+H)\^ 



2vr Jc 2 (A + iw)iw \ w(w — 2i(f)'(9i)) 

Equality (47) was obtained by making the change of variables — w = z and 
using that 7(^1) and 4>'(0i) are even and odd functions of 9\, respectively. 
In view of (47), in order to show that A{A) = — A{— A), it suffices to show 
that 



2tt Jc (A + iw)iw \ w(w — 2i<fi'(Qi)) 
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where C = C\ + C2 is the contour corresponding to the circle with radius 8. 
Now, 

If A f 2( 7 (fli) -j(9 1 +iw)) \ 



2tt Jc (A + iw)iw V w(w — 2i<fi'(9i)) 

1 ; 2vr 7_c - *A) V w(w-iA) J 

If A f w-iA \ 

(49) + 27 1c «,(«,- iA) l0g U-i20'(6» 1 )J ^ 

We will show that both terms (48) and (49) vanish. We first consider (49). 
For 7 € [0, 1] and a G [— <5, 6], define f{j) as 

1 f A 

/(7) = ^- / —, r-rz log(-fw - ia) dw. 

2tt J-c w(w — i A) 

Using residue calculus (see [14], page 224), it is easy to see that /(0) = 0. 
Indeed, a standard dominated convergence argument yields 



A 

dw = 0, 



C (w — iA)(jw — ia) 



where the previous integral has again been evaluated using residue calculus. 
As a result, we obtain that /(l) = 0. Applying these considerations with 
a = A and a = 2(f)' shows that the integral in (49) equals zero. We also 
can apply residue calculus to evaluate (48) directly as follows. Consider 

h (w)= , A lo /^i)-7(0i + H) 



w(w — iA) \ w(w — iA) 

Using the change of variables w = h + i A and the definition of A = 9\ — 
9$ with 7(6*1) = 7(00)1 we can evaluate the residue of /1 at w = iA as 
Residue(/i;iA) = — ilog(— 27' (#0)/ A). We also can obtain Residue(/i; 0) = 
ilog(27'(6 , i)/A). Therefore, using residue calculus, we obtain that the inte- 
gral in (48) equals 

-i log(-2 7 , (0 o )/(27 / (0i))) = -i log(7'(0i)/y (fli)) = 0, 

since in the case of symmetric distributions 7' (A) is odd and 9\ = —9q. 
Finally, we analyze B(A). Note that 

(50) B(A) = J- f " A log(2(l - g(\))\~ 2 ) d\ 

2tt J\\\>6 {A + t\)i\ 

(51) + V, J lM> _s {A + iXjiX bg I 1 " A-2#i)J ^ 
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Let Bi(A) and -82(A) be denned as (50) and (51), respectively. Since X\ is 
symmetric, it follows that log(2(l — g(X))X~ 2 ) is real. As a result, we obtain, 
just by integrating the real and imaginary parts of the integrand in B\(A), 



(52) Bi(a) = I L_ s ^ - S(A > )A "» 



Expression (52) yields an asymptotic expansion in odd powers of A for 
-Bi(A). Again, integrating the real and imaginary parts in 52(A), we obtain 



(53) ^'^LiA^^'^ 1 A-2 Wl ) 



The previous identity for 52(A) is obtained by observing that the integral 
of the imaginary part must vanish. This occurs because, for all 9\ small, 
the function log(l — iXH(9±, X)/(iX + 2<^'(#i))) satisfies the parity prop- 
erty, which can be verified by observing that, since 7(iA) = Eqcos(XX) + 
iEQsm(XX), it follows that hk(iX) £ V; also, using Proposition 6, we obtain 
that iX/(iX + 2<f>'{9{)) satisfies the parity property. Therefore, the closure 
properties proved in Proposition 6 together with an expansion of the loga- 
rithm yield that log(l - iXH(0 x ,X)/(iX + 20' (0i))) G V, which justifies (53) 
and (54). For notational convenience, let us define 



(55) C(9 1 ,X) = Y / h2k(iX)9 2 1 k /(2k)\ 

k=l 

and 

00 

(56) D(0 U A) = h 2k ^i(iX)9 2k - 1 /(2k - 1)!, 

k=l 

where hk(iX) = (j( k '(iX) — /Xfc)/(1 — j(iX)) — 2i^k+\/X. Since the distribution 
of X± is symmetric, we have that j(iX) is even and real. Moreover, we also 
have that hk(i\) is even if and only if k is even. We also can see that 
Re(H(9, A)) = C(9, A) and lm(#(0, A)) = D(9, A) are even and odd functions 
of both 9 and A [meaning that, for every 9 £ (—77/2, 77/2) fixed, C(9, •) is even 
and, similarly, for each A € R, C(-, A) is also even on (—77/2,77/2), say]. Using 
this notation, we can write 

XH(9i,X) X 2 C{9 l ,X) - 2(f)'(9i)XD(9i,X) 



(57) 



X-2<f>'(9 1 )i X 2 + {24>>{9 1 )) 2 

.2(f>'(9 l )XC{9 1 ,X) + X 2 D{9 1 ,X) 



(58) +i- 



X 2 + (2<P'(9 1 )Y 
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Let us define C(9\,X) and D(9i,X) as the real and imaginary parts of 
XH(6i, A)/(A — 2(f)' (6i)i), respectively, as indicated in the corresponding ex- 
pressions (57) and (58). Since XH(6i, A)/(A — 2(p'(9i)i) holds the parity prop- 
erty, C(9\, A) and D(6\, A) are even and odd function in both arguments 9\ 
and A. By symmetry of the distribution of X±, we have that A = 29\, also as 
a consequence of symmetry, 2<f)' (9±) is an odd (real) analytic function of 9\ 
at the origin, which implies that (2(//(#i)) 2 is even. Hence, using the expan- 
sion of log(l — z) at z = in expressions (53) and (54) (justified by virtue 
of Proposition 3), we see that an asymptotic expansion for the integral (53) 
involves expanding expressions of the form 

(59) K(0i) =~ / TX ^r C(9 1 ,X) k D(e 1 ,X) 2m dX, 

2ir J\\\>s (A z + X z ) 

where K(9±) is an even function of 9\ which is also (real) analytic at the 
origin. This implies, in view of (55) to (58) and the properties of 2<//(# 1 ) 
discussed before, that an asymptotic expansion for (59) must be given in 
odd powers of A only, which must be also the case for the integral in (53). 
The treatment for the integral (54) is completely analogous and also yields 
an asymptotic expansion in odd powers of A. This yields the conclusion of 
the theorem. □ 
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